Accuracy of unperturbed motion of particles in a gyrokinetic semi-Lagrangian code
نویسندگان
چکیده
Inaccurate description of the equilibrium can yield to spurious effects in gyrokinetic turbulence simulations. Also, the Vlasov solver and time integration schemes impact the conservation of physical quantities, especially in long-term simulations. Equilibrium and Vlasov solver have to be tuned in order to preserve constant states (equilibrium) and to provide good conservation property along time (mass to begin with). Several illustrative simple test cases are given to show typical spurious effects that one can observes for poor settings. We explain why Forward Semi-Lagrangian scheme bring us some benefits. Some toroidal and cylindrical GYSELA runs are shown that use FSL. Key-words: gyrokinetic, conservation laws, FSL ∗ CEA Cadarache, IRFM bat. 513, 13108 Saint-Paul-les-Durance † University of Strasbourg, 7 rue Descartes, 67000 Strasbourg & INRIA CALVI ‡ INRIA IPSO & IRMAR, 263 av Général Leclerc, 35042 Rennes ha l-0 07 27 11 8, v er si on 1 3 Se p 20 12 Précision du mouvement de particules non perturbées dans un code gyrocinétique Résumé : Une description imparfaite de l’équilibre peut conduire à des artéfacts numériques dans des simulations turbulentes gyrocinétiques. Aussi, le solveur de Vlasov and le schéma d’intégration en temps ont un impact fort sur la conservation de quantités physiques, notamment lorsque les simulations sont en temps long. L’équilibre et le solveur Vlasov doivent être finement choisis et paramétrés pour conserver les états constants (équilibre) et pour autoriser de bonnes propriétés de conservation en temps (en commençant par la conservation de la masse). Plusieurs cas tests illustratifs sont donnés pour montrer les problèmes numériques typiques que l’on peut observer si les choix pris ne sont pas adéquats. Nous expliquons pourquoi le schéma FSL (Forward Semi-Lagrangian) apporte une réponse à certains problèmes. Des simulations en configuration cylindrique et torique de GYSELA sont présentées qui utilisent FSL. Mots-clés : gyrocinétique, lois de conservation, FSL ha l-0 07 27 11 8, v er si on 1 3 Se p 20 12 Unperturbed motion of particles 3
منابع مشابه
A Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations∗)
A semi-Lagrangian scheme is applied for the first time to computations of charged particle motions along magnetic field lines, to numerically solve the δf gyrokinetic equations in a flux tube geometry. This new solver adopted in the gyrokinetic Vlasov simulations has an advantage over the conventional Eulerian codes in calculating the parallel dynamics, because semi-Lagrangian schemes are free ...
متن کاملGYSELA , a full- f global gyrokinetic Semi-Lagrangian code for ITG turbulence simulations
This work addresses non-linear global gyrokinetic simulations of ion temperature gradient (ITG) driven turbulence with the GYSELA code. The particularity of GYSELA code is to use a fixed grid with a Semi-Lagrangian (SL) scheme and this for the entire distribution function. The 4D non-linear drift-kinetic version of the code already showns the interest of such a SL method which exhibits good pro...
متن کاملLIGKA: A linear gyrokinetic code for the description of background kinetic and fast particle effects on the MHD stability in tokamaks
In a plasma with a population of super-thermal particles generated by heating or fusion processes, kinetic effects can lead to the additional destabilisation of MHD modes or even to additional energetic particle modes. In order to describe these modes, a new linear gyrokinetic MHD code has been developed and tested, LIGKA (Linear Gyrokinetic Shear Alfvén Physics) [1; 2], based on a gyrokinetic ...
متن کاملA 5D gyrokinetic full-f global semi-Lagrangian code for flux-driven ion turbulence simulations
This paper addresses non-linear gyrokinetic simulations of ion temperature gradient (ITG) turbulence in tokamak plasmas. The electrostatic Gysela code is one of the few international 5D gyrokinetic codes able to perform global, full-f and flux-driven simulations. Its has also the numerical originality of being based on a semi-Lagrangian (SL) method. This reference paper for the Gysela code pres...
متن کاملImproving conservation properties of a 5D gyrokinetic semi-Lagrangian code
In gyrokinetic turbulent simulations, the knowledge of some stationary states can help reducing numerical artifacts. Considering long-term simulations, the qualities of the Vlasov solver and of the radial boundary conditions have an impact on the conservation properties. In order to improve mass and energy conservation mainly, the following methods are investigated: fix the radial boundary cond...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- CoRR
دوره abs/1209.0317 شماره
صفحات -
تاریخ انتشار 2012